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Nonlinear dynamics of wave packets in T^T-symmetric optical lattices near the phase-transition 
point are analytically studied. A nonlinear Klein- Gordon equation is derived for the envelope of 
these wave packets. A variety of novel phenomena known to exist in this envelope equation are 
shown to also exist in the full equation including wave blowup, periodic bound states and solitary 
wave solutions. 



Linear Schrodinger operators with complex but parity- 
time (7^T)-symmetric potentials have the unintuitive 
property that their spectra can be completely real yp. 
This phenomenon was first studied in quantum mechan- 
ics where a real spectrum is required to guarantee real 
energy levels. The same phenomenon was later investi- 
gated in paraxial optics, where T^T-symmetric potentials 
could be realized by employing symmetric index guid- 
ing and an antisymmetric gain/loss profile [2l|3]. In this 
optical setting, VT potentials have been experimentally 
realized [4, 5j. In temporal optics, T^T-symmetric lat- 
tices have been experimentally obtained as well [6]. So 
far, a number of novel physical phenomena in optical 
VT systems have been reported, including phase tran- 
sition (T^T-symmetry breaking), nonreciprocal Bloch os- 
cillation, unidirectional propagation, distinct pattern of 
diffraction, formation of soliton families, and so on 0]- 
[TT] . The search of additional new behaviors in optical 
VT systems is still ongoing. 

In this Letter, we analytically study nonlinear dynam- 
ics of wave packets in T^T-symmetric optical lattices near 
the phase-transition point (where bandgaps close and 
Bloch bands intersect transversely like the letter 'x'). A 
nonlinear Klein-Gordon equation is derived for the en- 
velope of these wave packets near the band intersection. 
Based on this envelope equation, we predict a variety of 
novel phenomena such as wave splitting, wave blowup, 
periodic bound states and solitary wave states. We fur- 
ther show these predicted phenomena occur in the full 
model as well. 

The paraxial model for nonlinear propagation of light 
beams in T^T-symmetric optical lattices is 

i^. + + V{x)^ + = 0, (1) 

where z is the propagation axis, x is the transverse axis, 

V{x) = [cos(2x) + iWo sin(2x)] (2) 

is a T^T-symmetric potential, Vq is the potential depth, 
1^0 is the relative gain/loss strength, and a = ±1 
is the sign of nonlinearity. All variables are non- 
dimensionalized. 

First we discuss the linear diffraction relation of Eq. 
([T]) at the phase-transition point Wq = 10 E] • In this 
case, the linear equation of (fTl) can be solved exactly pT] . 



Its Bloch solutions are 

*±(x,z;m) = I±(k+2m) (V^oe*") e-'^\ (3) 

where Ik is the modified Bessel function, /i = (/c + 2m) ^ 
is the diffraction relation, k is in the first Brillouin zone 
k e [—1,1], and m is any nonnegative integer. This 
diffraction relation is depicted in Fig. [ijA), where dif- 
ferent colors indicate different Bloch bands. At /c = 
and ±1, adjacent Bloch bands intersect each other trans- 
versely like the letter 'x'. Due to this intersection (degen- 
eracy), wave packets near these points will exhibit novel 
dynamics which we will reveal next. 

At these intersection points, ^+ = thus Bloch 
solutions ([3| are degenerate and 27r-periodic in x. Posed 
as an eigenvalue problem for ^ = (j){x)e~^^^ in the linear 
equation of ([T]), we get Lcj) = — where L = Oxx-^Voix), 
and Vo{x) is the VT lattice at the phase-transition point 
Wo = 1. Then at these band-intersection points, the 
eigenvalues are ji = n? ^ where n is any positive integer. 
The corresponding eigenfunctions are 

^M^7..(Pie-) = f (4) 

These eigenvalues ji = n? all have geometric multiplicity 
1 and algebraic multiplicity 2, thus there exist a general- 
ized eigenfunction cj)^ satisfying 

{L + ii)4>' = 4>. (5) 

For Vb = >/6 and n = 1, this eigenfunction and general- 
ized eigenfunction are plotted in Fig.[ljB). 

Now we consider nonlinear dynamics of wave pack- 
ets near these band-intersection points. For convenience, 
we first take the VT lattice to be exactly at the phase- 
transition point (i.e., = 1). Generalization to lattices 
near the phase-transition point will be made afterwards. 

Nonlinear wave-packet solutions near a band- 
intersection point can be expanded into a perturbation 
series 

^ = e-'^'[eA{X, Z)^{x) + e^^i + e^V^2 + -..], (6) 

where ji = v? is the propagation constant at the band 
intersection, (/)(x) is the degenerate Bloch mode given 
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(A) Dispersion Relation 



(B) Degenerate Bloch Mode 



where 




FIG. 1: (A) Diffraction relation; (B) Magnitude and phase of 
eigenfunction (solid blue) and generalized eigenfunction 0^ 
(dashed red) for fi — 1 and Vb = V^; (C) Linear unidirectional 
wavepacket and (D) linear wavepacket splitting in Eq. ([T]) at 
the phase-transition point Wo = 1. 



in Eq. Q, A{X^Z) is the envelope of this Bloch mode, 
X = ex^ Z = ez are slow variables, and < e <C 1 
is a small positive parameter. Substituting the above 
perturbation series into the original equation ([T]), this 
equation at 0(e) is automatically satisfied. At order 
we have 

(L + /i)V^i = -iAz^ - 2Ax^:.. (7) 

The solution to this equation is 

V'l = -iAz<l>' - 2Ax<t>\ (8) 

where 0^ is the generalized Bloch mode defined in (j5]), 
and 



(9) 



This (j)^ solution exists and can be determined by Fourier 
series [11 . 

At O(e^) we get 

(L + /i)V^2 = -Azz^' + Axx - 0) 

Vi2Azx (0' + - cj\A\^A\c^\^. (10) 

The solvability condition of this equation is that its right 
hand side be orthogonal to the adjoint homogeneous so- 
lution 0*. Using Fourier expansions of solutions (j)^ and 
(j)^ (see [11] ) , this solvability condition then yields the fol- 
lowing nonlinear Klein-Gordon equation for the envelope 
function A{X, Z): 



7 = (-l) 



n+l 



\(t)\^(t)^dx. 



(12) 



In view of the formula Q for (/), it is easy to see that 
sgn(7) = {—l)^~^^a. For the values of Vq = V6 and 
n = 1 which we will use in later numerical simulations, 
7 ^ 11.0430cr. 

Now we extend the above envelope equation to the case 
where the VT lattice is near the phase-transition point 
(i.e.. Wo ~ 1). Following similar perturbation analy- 
sis, we find that when ja = n = 1 (the lowest band- 
intersection point) and Wq = 1 — ce^, the envelope is gov- 
erned by a slightly more general nonlinear Klein- Gordon 
equation. 



Azz - ^ri^Axx + + -i\AfA = 0, 



(13) 



where 7 is as given in (12), a = and the i/ji solu- 

tion in ^ is still given by Eq. ([8|. When 11 = 11? with 



n > 1 and Wq = 1 — ce, the envelope equation ( |13| ) will 
contain an additional term proportional to lAz- But this 
lAz term can be eliminated through a gauge transforma- 
tion A Ae~^^^o -1)^ aj^d the transformed equation 
remains the same as (13), except that a = c^V§ /QA for 
n = 2 and a = for n > 2. If c = (i.e., at the phase- 
transition point), then a = 0, hence Eq. (IT3| reproduces 



Eq. (11) as a special case. When n > 1, the solution 
([8| will also contain an additional term proportional to 
icA. 



The nonlinear Klein- Gordon equation (13) is second- 
order in Z, which means that two initial conditions, 
A(X, 0) and Az{X, 0), are needed. These two initial con- 
ditions can be obtained from the initial conditions of the 
perturbation series ^ at orders e and e^, i.e., from the 
initial envelope A{X,0) as well as V^i|2=o- This V^i|2=o 
generally contains many eigenmodes of the operator L, 
but only the (j)^ component in it affects the dynamics of 
envelope A. If we denote B{X) as the envelope function 
of the (j)^ component in V^i|2=o, then by projecting tlJi\z=o 
(such as from ([8|) onto (/)^, we find that 



B{X) = -i[Az{X,0)+2nAx{X,0)]. (14) 

This relation holds for all n values. Thus, initial con- 
ditions for the original VT model ([T]) and those for the 
envelope equation (13) can be related as 



^(x, 0) = eA{X, 0)(^(x) + e^B{X)(l)^{x). (15) 

This initial- value connection will be used in our numeri- 
cal simulations later. 

Now we examine the envelope dynamics in the non- 
linear Klein-Gordon equation (13), and show that the 



Azz-^n^Axx^7\A\^A = 0, 



(11) 



corresponding wavepacket dynamics occurs in the origi- 
nal VT model ([T]) too. In this discussion, we take n = 1, 
Vb = >/6 and e = 0.1 for definiteness. 

First we consider the linear Klein-Gordon equation 
(13) at the phase-transition point, i.e., Azz —^n^Axx = 
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0. This is the famihar second-order wave equation. It ad- 
mits unidirectional wave solutions F{X — 2nZ) as well as 
bidirectional wave solutions F{X — 2nZ) + G{X + 2nZ). 
In the original linear VT model ([T]), we find that the 
corresponding wavepacket solutions also exist. Examples 
are displayed in Fig.[TJC,D). Note that the wave splitting 
in Fig.jljD) is manifestation of the 'x'-shaped diffraction 
structure at the band-intersection point, while the unidi- 
rectional propagation in Fig. [TJC) occurs for B{X) = 0. 



Envelope \A\ Full solution |^| 




X X 



FIG. 2: Nonlinear wavepacket solutions below the phase- 
transition point under self-defocusing nonlinearity. Top: a 
blowup solution; bottom: a periodic bound state. Left: enve- 
lope solutions \A\; right: full solutions |^|. 

Next we consider envelope solutions in the nonlinear 
Klein-Gordon equation ([l3| below the phase-transition 
point, i.e., Wq < 1, or a > (above the phase-transition 
point, infinitesimal linear waves are unstable, thus it is 
not pursued). When the nonlinearity is self-defocusing 
(cr = -1), 7 ^ -11.0430 < 0. In this case, en- 
velope solutions can blow up to infinity in finite dis- 
tance pj]. One such example is shown in Fig. |2] (up- 



per left panel) for c = 1 (a = 18) and initial condi- 
tions A(X,0) = 1.2sech(X), Az{X,0) = -2Ax(X,0) 
{B{X) = 0). In the fuh nonlinear VT model Q, we 
have found similar blowup solutions which are displayed 
in Fig. [2] (upper right panel). This solution-blowup un- 
der self-defocusing nonlinearity is very surprising. Note 
that in the full model ([T]), our asymptotic envelope ap- 
proximation breaks down as the singular (blowup) point 
is approached. In this case, the amplitude of the full- 
model solution eventually saturates, but its power still 
grows unbounded pT| . 

Under self-defocusing nonlinearity, the envelope equa- 
tion ([l3| also admits X-localized and Z-periodic bound 
states. One example with c = 1 (below the phase- 
transition point) is shown in Fig. [2] (lower left panel). The 
initial condition for this solution is A(X, 0) = sech(X) 
and Az{X,0) = {B{X) = -2iAx{X,0)). In the full 
nonlinear VT model ([T]), we have found similar periodic 
bound states as well, see Fig. [2] (lower right panel). 

Under self-defocusing nonlinearity and below the 
phase-transition point, the envelope equation ([T3| also 
admits stationary solitary waves A(X, Z) = F(X)e^^^ 
when uj lies inside the bandgap —^/a < u < ^/a. These 
solitary envelope solutions correspond to the solitons of 
the full nonlinear VT model ([T]) reported in [71 [11] . 

If the nonlinearity is self- focusing {a = 1), envelope 
solutions do not blow up, periodic bound states cannot 
be found, and stationary solitary waves do not exist in 
the envelope equation. In this case, nonlinear diffracting 
solutions similar to the linear diffracting pattern reported 
in [9 exist. In addition, solutions periodic in both X and 
Z can be found. 

In summary, we have shown that the nonlinear Klein- 
Gordon equation governs the envelope dynamics of 
wavepackets in VT lattices near the phase-transition 
point. We have also shown that a variety of novel phe- 
nomena in this envelope equation (such as wave blowup 
and periodic bound states) occur in the full VT model 
too. These findings open new possibilities for wave engi- 
neering in VT lattices. 
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